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Abstract 

We study the twist-three fragmentation function contribution to the single transverse spin asym- 
metries in inclusive hadron production in pp collisions, p^p — > h + X. In particular, we evaluate 
the so-called derivative contribution which dominates the spin asymmetry in the forward direction 
of the polarized proton. With certain parametrizations for the twist-three fragmentation function, 
we estimate its contribution to the asymmetry of 7r° production at RHIC energy. We find that the 
contribution is sizable and might be responsible for the big difference between the asymmetries in 
n and 7r° productions observed by the STAR collaboration at RHIC. 
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Single-transverse spin asymmetries (SSAs) in hadronic processes, such as in the sin- 
gle inclusive hadron production in single transversely-polarized nucleon-nucleon scatter- 
ing, p^p — > hX, have attracted much interests from both experimental and theoretical 



sides in the last 
derlying physics 



ew years, and great progress has been made in understanding the un- 
it]. Although it is a simple observable, defined as the spin asymme- 
try when one flips the transverse spin of one of the hadrons involved in the scattering: 
Ajsi = (da(S±) — da(— S±))/ (da(S±) + dcr(—S±)), it is far more complicated to explain in 
the fundamental theory of strong interaction. It usually represents a correlation between 
the transverse polarization vector S± of one of the hadrons and the transverse momentum 
Ph± °f the final-state hadron in the differential cross section. For example, in the process 
p^p hX, it is the correlation between the polarization vector S± of the incoming nucleon 
and the transverse momentum Ph± of the final-state hadron that generates the asymmetry. 
In this paper, we will focus on the SSAs in these processes, especially for the neutral mesons, 
it and 7] production, motivated by the recent striking experimental observations of large 
SSAs for them by the STAR collaboration at RHIC experiments [2, [3]. In particular, at 
forward direction of the polarized proton, it was found that the SSA for r] meson is much 
larger than that for 7r° meson. This also confirmed previous experimental observations from 
the fixed target experiments \&\. In addition to these STAR results, both BRAHMS and 
PHENIX collaborations at RHIC have observed large single transverse spin asymmetries in 
charged and neutral meson production in the forward rapidity region {5, 6]. 

In the QCD framework, there have been mainly two approaches to study the SSAs in high 
energy scattering processes: the transverse momentum dependent (TMD) parton distribu- 



tion and fragmentation function approach 



approach in the collinear factorization 
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15] and the twist-three quark-gluon correlation 



1 6M2 lj . In the TMD approach, it is required to have 



an additional hard momentum scale besides the transverse momentum Ph± of the hadron, 
such as Q 2 , the momentum transfer square of the virtual photon in semi-inclusive deep in- 
elastic scattering (SIDIS) or Drell-Yan lepton pair production in pp collisions, to study the 
transverse momentum dependence in these processes. On the other hand, the twist-three 
approach is more appropriate in the processes that all momentum scales are much larger 
than the nonperturbative scale Aqcd [12|- It has been shown that these two approaches 
are consistent with each other for the SIDIS and Drell-Yan processes in the intermediate 



transverse momentum region where they both apply 22 
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The TMD approach has also been used to calculate the SSA in the single inclusive hadron 



production in p^p collisions in a model dependent way 10|, |24], |25|. However, there has been 
no factorization argument based on TMD parton distribution for this process 26|. In this 
paper, we follow the collinear factorization approach to study the single spin asymmetry 
coming from the twist-three fragmentation function contribution. In this process, there is 
only one large momentum scale, the transverse momentum Ph± of the final-state hadron, 
and the SSA is naturally suppressed by 1/Phx at large Ph± as a higher-twist effect, although 
how large this behavior should manifest is not clear. At low transverse momentum, it 
will ultimately fail, and a modification has to be made to account for the experimental 
data 

In the twist-three collinear factorization approach, for the SSAs in hadron production in 

p\A)+p{B)^h + X , (1) 

the twist-three effects come from the distribution functions of the incoming polarized nucleon 
(A) with momentum Pa or the unpolarized nucleon (B) with momentum Pb, or the fragmen- 
tation function for the final state hadron (h) with momentum P^. Therefore, schematically, 
we can write down the single transverse spin dependent differential cross section in the 
following way 

HQ, 



da(S±) = e < fS± a P h p J [dx][dy][dz] {<pf/ A (x, x') <g> <f)j/B{y) ® D h/c (z) <g> Uif\ c (x, x', y, z) 
+<f> i/A {x) ® <f>f} B (y, y') <g> D h/c {z) ® U^Jx, y, y', z) 

+4> i/A (x) <g> <f> j/B (y) ® D$ e (z, z') ® H^ c (x, y, z, z') ) , (2) 



where is defined as e*f = e^ va ^ P^Pbv / Pa ■ Pb with convention of e 0123 = 1. Here we 
use _L to label the transverse direction in the center of mass frame of P A and Pb- In the 
above equation, ® stands for the convolution in the longitudinal momentum fractions x, y, 
and z. The leading twist parton distributions are labeled by 4>i/ A (x) and 4>j/B(y)- Because 
hadron A is transversely polarized, we immediately see that <pi/ A represents the leading-twist 
quark transversity distributions. However, for unpolarized hadron B, <f>j/B represents both 
leading twist quark and gluon distributions. Leading twist parton fragmentation function is 
represented by D h / C (z) where parton c can be a quark or gluon. In the above equation, the 
superscript (3) represents the twist-three correlations for the distribution functions or the 
fragmentation functions. Since these twist-three functions normally involve two variables in 



longitudinal momentum fractions, we have made them explicit in the above formula. For a 
complete analysis, the above three terms have to be taken into account. The first term in the 



above equation, the contribution from the twist-three 



parton distributions 



x, x') from 



the polarized nucleon have been calculated in Refs. 17H19| . The second term from the twist- 
three effect in the unpolarized nucleon cj)^ 3 ) B (x,x') is found very small in the forward region 



of the polarized nucleon where the largest asymmetry has been found experimentally 28]. 
The third term is least known in the literature, which we will focus in this paper. There have 
been attempts to formulate this part of contribution 29j. However, a universality argument 



for the Collins fragmentation function 



30 



lated in Ref. 



would indicate that the contribution calcu- 



29] vanishes. This universality has also been extended to pp collisions 32 ]. 



Large azimuthal asymmetries coming from this Collins effect have been observed in semi- 
inclusive hadron production in deep inelastic scattering and di-hadron production in e + e~ 
annihilation processes by the HERMES [33J] and COMPASS [34j, and BELLE collabora- 
tions 35)], respectively. In a recent publication, two of us have re-examined the twist-three 
fragmentation contribution to the single spin asymmetry [23] . In particular, we have identi- 
fied the twist-three fragmentation function corresponding to the TMD Collins fragmentation 
function, and shown that the TMD and collinear factorization approaches are consistent in 
the intermediate transverse momentum region in the SIDIS process. In this paper, we will 
extend this formalism to the twist-three fragmentation function contribution to the above 
process of ([I]), the single spin asymmetries in inclusive hadron production in p^p scattering. 

The twist-three fragmentation function which corresponds to the universal Collins TMD 
fragmentation function, is defined as 1 



H(z) = — 
v ; 2 



2tt 2 
x (P h X\ij(0)\0) + h.c.} 



'D'r-r / d( ///••" 'k 



m\p h x) 



(3) 



where k + = P^/z. This function can also be written as the transverse momentum moment 
of the Collins fragmentation function H^-(z,p^), 

\Pt\ 2 



Hz) 



(4) 



In the above definition Eq. 03]), we have set the hadron momentum in the z direction. The 
index a is in the transverse direction perpendicular to the hadron momentum direction, and 



1 In the definition of H (z) in Ref. 23j, a factor of 1/2 missed 
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is not considered to be summed up. In order to uniquely define this transverse component 
(index), we introduce a light-like vector oc (P°, —Ph) with normalization that rih-Ph = l 2 - 
Here and in the following, we use subscript "T" to distinguish the transverse direction 
perpendicular to Ph from "_L" in Eq. (2) for transverse direction perpendicular to Pa and 
Pb- Therefore, the transverse component for any vector p^ can be defined as pt^ = p^ — p ■ 
rihPh — p ■ Phfih- F rom this definition, we can immediately see that px is space-like, and 
perpendicular to P^: pr ■ Ph = 0. 

As discussed in Ref. [23], the above twist-three fragmentation function belongs to more 
general two-variable dependent twist-three fragmentation function defined as, 

M*,**) = Z ~fJ ^^e^e^<-\{Tra^(0\im(Cmr)\^ 

x (P h X$(0)\0) + h.c.} , (5) 

where kf = P^ / z i an d kg = — k 2 . Similarly, we can define a F-type fragmentation 
function H F (zi, z 2 ) by replacing D% with F +a in Eq. (J5j). -D-type and F-type functions are 
related to each other, 

H D ( Zl ,z 2 ) = Pvl^— r )H F (z 1 ,z 2 ) + 5(---)^H(z 1 ). (6) 

From this equation, we can see that Ho is more singular than Hp at Z\ = z 2 - Moreover, 
recent study has shown that Hp(zi,z 2 ) vanishes when z\ = z 2 [31] . Therefore, it is more 
convenient to express the cross section in terms of H(z) and Hp(zi, z 2 ), which we will follow. 

In this paper, we are interested in obtaining the phenomenological important contribu- 
tions in hadron production in the forward region of the polarized beam. Similar to the 
twist-three distribution function contributions, we find that the twist-three fragmentation 
function contributions also contain a derivative term which certainly will dominate this part 
of contribution in the forward rapidity region. In addition, the derivative contribution from 
the twist-three fragmentation function is associated with H(z) defined in Eq. ([3]), whereas 
the contribution associated with H F (zi, z 2 ) does not contain any derivative contributions. In 

2 In the low transverse momentum semi-inclusive hadron production in DIS process, will be proportional 
to the incoming nucleon momentum as we used in Ref. [23] . In current study, there is no natural available 
momentum for this choice. We can also perform the calculations without choosing the vector n^, but 
keeping all transverse component in Eq. ([3]) perpendicular to the hadron momentum. This will lead to 
the same results. 
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FIG. 1: Generic Feynman diagrams to calculate the derivative term of the twist-three fragmentation 
function contributions to the single inclusive hadron production inpp scattering, p^p — > hX: p\ and 
P2 are the two incident partons ' momenta, k\ and ki are the outgoing partons ' momenta, where the 
quark k\ fragments into final-state hadron Ph- The expansion of the scattering amplitude in terms 
of transverse momentum component of ki = Ph/z + k\T leads to the twist-three contributions from 
the fragmentation function H{z) defined in Eq. 

this paper, we will report the derivative-term results, and carry out the numerical estimates 
for their contributions to the single spin asymmetries at RHIC. We will leave the detailed 
derivations, including the derivative and non-derivative terms in a future publication. 



We follow the technique developed in the last few years [l7H23l| to derive the contributions 
from the twist-three fragmentation functions. First, we calculate the associated scattering 
amplitudes in terms of various twist-three matrix elements of hadrons. From diagrammatic 
point of view, these individual contributions are not in a gauge invariant way. However, 
the final results will be gauge invariant when we sum up all the contributions. In terms 
of twist-three operators, we have the contributions from the Ot^P-, At, and 8tA + (d + At)- 
These contributions indeed will form the gauge invariant results in terms of H(z), Hf as 
defined above. In particular, drip term corresponds to H(z) whereas At term corresponds 
to Hp. Since the derivative terms are associated with the function H(z), we will need to 
calculate the contributions from the Dt^ part. In order to carry out this calculation, we have 
to perform the collinear expansion for the hard partonic part associated with this matrix 
element. In Fig.IU, we illustrate the generic Feynman diagrams contributing to the derivative 
terms. In this diagram, pi and P2 are the two incoming partons' momenta from the polarized 
and unpolarized nucleons, respectively; k\ and k<i are the two outgoing partons' momenta 
where k\ fragments into the final-state hadron Ph- The derivative terms come from the 



6 



transverse momentum expansion of hard partonic scattering amplitudes, in particular, from 
the on-shell condition for the unobserved particle k 2 in Fig. [TJ This momentum can be 
written as k 2 = p\ + p 2 — ki, where we can parameterize k\ as k\ = Ph/z + k\j>. Thus, the 
on-shell condition will lead to the following expansion, 

S((pi +P2~ hf) = S{(px + p 2 - P h /zf) - 2(pi + pa) • k 1T 5'(( Pl +p 2 - P h /zf) , (7) 

where the second term will result into a derivative term contribution to the transverse 
spin dependent differential cross section. As discussed above, for convenience, we use 
a light-like vector to constrain the transverse momentum expansion which is per- 
pendicular to the final state hadron's momentum. Therefore, kyr can be expressed as 
k^ T = ki — k\ ■ PbJ^ h — ki • rihP^. Substituting this expression into the above equation, 
we will obtain the contribution from the matrix element associated with the operator dxtp 
which will result into the contribution of H(z). The derivative on the Delta function will 
be translated into the derivative on the twist-three fragmentation function H(z) by partial 
integral. This contribution is similar to what has been calculated for the twist-three distri- 
bution functions 17| . In particular, the hard part is proportional to that from the relevant 
twist-two Born diagrams. In our case, it is the transverse spin transfer coefficients calculated 
in Ref. 36]. However, there exist major difference between these two contributions. The 
twist-three distribution contributions come from the pole contributions from the initial/final 
state interactions, and the derivative terms have contributions from the expansion of the 
on-shell condition for the unobserved particle and the double pole contribution from the final 
state interaction diagrams [17|. In the twist-three fragmentation function contributions, in 
contrast, we do not take pole contributions, and the derivative terms only come from the 
expansion of the on-shell condition of the unobserved particle in the above 2 — > 2 processes. 
The final result for the derivative contributions of H (z) takes the following form, 

d ( H(z) 



d 3 Aa(s ± ) Qa 2a 2 sX ^ f 1 dx> 1 f 1 dz 

- «• C X min J Z mtn 

P£ + z( P2 .P h ni-p 2 .n h P£ 



where 




-z- 



—zu 



dz \ z 2 

H ab ^ c (s,i, u) 



S + U/z' mm S 



T + U 

1 j 



(9) 
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and the hard factors H a h^ r are defined as 



H. 



qq-^qq 



Hqq'—^qqi 



AN 2 



N 2 



1 Asu 



H, 



qq-^qq 



qq^tqq 



H, 



Hqq^qq 
Hqq^qg 

N 2 



AN 2 -f 2 

Asu 1 4s 
Asu 



-i 2 

N 2 -l 



AN? 



1 Au 

r 2 + iVcT 



N 2 



qg^qg 



1 1 Asu 



(10) 



where s, t, and u are the usual partonic Mandelstam variables. They are the same as the 



transverse spin transfer hard coefficients calculated in Ref. 



36(. Substituting the expression 



of rih into Eq. (8), the factor in the bracket can also be written as, 



— >• 



-zu 



z I x(-u) + x'(-t) 



(11) 



In the forward rapidity region of the polarized nucleon, we have x x' and — u —t, and 
we can further simplify the transverse spin dependent differential cross section as 



_ d 3 Aa(s ± ). p 2a 2 s ^ f 1 dx' , 1 f 1 dz 

— cFp — \f° rward ~ Mi^yL /, ~j^ x >~ ha y x > / ~ 

a.b.C X 'min J Zmin 



d ( H{z) 
dz 



x- 



-H ab ^Js,t,u) . 



(12) 



x'S + T/z —zu 

This term will be most phenomenological relevant for the single spin asymmetries of hadron 
production in the forward direction of the polarized nucleon. 

In order to demonstrate the twist-three fragmentation function contribution to the SSA 
in inclusive hadron production in p^p collisions, we need the unknown, but universal, twist- 
three fragmentation function H(z). We notice that H(z) can be related to the Collins 
fragmentation function H^z^p 2 ^) as in Eq. which has been studied from the available 



experimental data 



37, 



]. However, we emphasize that the Collins function H^(z,p\) are 
fitted from small transverse momentum region where TMD factorization applies. To obtain 
the functional form for H(z), one has to assume a transverse momentum dependence in the 
Collins fragmentation functions. In principle, the twist-three fragmentation functions H(z) 
should be extracted from the experimental data of the SSAs at large transverse momentum 
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region where collinear factorization applies, similar to what has been done in [18j, or from 
the transverse momentum weighted azimuthal asymmetry measurements in which the twist- 
three fragmentation function H{z) enters directly. 

For the purpose of estimating its contribution to the SSAs in the inclusive ir° production 
in p^p collisions at RHIC energy, we parameterize the twist-three fragmentation function in 
terms of the leading-twist fragmentation function, 

H{z) = C f z a D(z) . (13) 

We leave a more comprehensive parameterization for these fragmentation functions in a 
future study, including that for charged mesons as well. The additional z factor in the 
parameterization comes from the consideration that this novel fragmentation is mostly a 
valence behavior. Normally, we would also have a (1 — z) suppression for the twist-three 
function. However, for scalar hadron production, the twist-three fragmentation function 
is not power suppressed in terms of (1 — z), similar to the power counting of the Boer- 



Mulders function of 7r meson at large- a; 



39J. We have also adopted the quark transversity 



distributions from the parameterizations in Ref. |40J and the unpolarized fragmentation 
function in [41]. In Fig. [21 we show the predictions with the above parameterizations for the 
it production with Cf — —0.4. The three curves from up to bottom correspond to a = 1 
(solid), a = 2 (dashed), and a = 4 (dotted), respectively. With our parametrization of H(z), 
the twist-three fragmentation function can generate a sizable SSA in inclusive tt° production 
at RHIC energy y/s = 200 GeV. These contribution are comparable to that of the twist- 



three distribution from the polarized nucleon which have similar parameterizations 



The single transverse spin asymmetry of r] meson has also been studied by the STAR 
collaboration at RHIC recently [2]. A significantly larger asymmetry An has been observed 
for i] meson compared to 7r°. As we discussed, in the twist-three collinear factorization, the 
dominate contribution comes from the twist-three distribution of the polarized nucleon, as 
shown in the first term of Eq. (J2j), and the twist-three fragmentation function, as shown in 
the third term of Eq. fl2]) . The difference between 77 and 7r° from the first term will be mainly 
due to the difference in the unpolarized meson fragmentation function. However, from the 
measurements of spin-averaged cross section of rj and 7r° by PHENIX collaboration at RHIC 



42J, one finds that for a large range of transverse momentum, tj/tt ratio is a constant. This 
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1.5 2 2.5 3 3.5 4 4.5 5 5.5 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

P hl X F 

FIG. 2: Twist-three fragmentation contributions to the single spin asymmetries in hadron produc- 
tion in the forward direction of the polarized nucleon at RHIC at y/s = 200GeV , as functions of 
xp, with the parameterization for the twist-three fragmentation function H(z) in Eq. i!3\) with 
parameters Cf = —0.4 GeV and a = 1,2,4 (from up to bottom), respectively. 

indicates that the unpolarized fragmentation function for rj and tt° is very similar, and will 
lead to the similar SSAs for them if this contribution dominates. In other words, it will be 
very difficult to explain the large difference between the SSAs of r] and 7r° mesons from the 
twist-three distribution contribution from the polarized nucleon 3 . Since the second term in 
Eq. generates a very small asymmetry, one would expect the large difference between r\ 
and 7T° would come from the twist-three fragmentation function contribution if one believes 
the collinear factorization is indeed the right approach to describe the SSA of inclusive 
hadron production in pp collisions. Since the twist-three fragmentation function H(z) for 77 
and 7T° in general need not to be the same, this might generate the needed difference observed 
by the experiments if r\ meson has a much larger twist-three fragmentation function H(z) 
compared to 7r°. To test this scenario, it will be very important to study the associated 
Collins fragmentation for 77 meson in e + e~ annihilation and/or semi- inclusive DIS processes 
and compare to that for 7r° meson. We hope that, in particular, the BELLE collaboration 
can carry out this measurement and cross check with the STAR observation. Meanwhile, 
we emphasize that to better understand the single spin asymmetry for these processes and 
finally pin down the difference, one needs to take into account both twist-three contributions 

It might still be possible that the strange quark contribution from the polarized nucleon may dominates 
and leads to a larger SSA in n meson production, which is, however, unlikely in the forward rapidity 
region (the valence region) of the polarized nucleon. 
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and perform a global analysis of the experimental data. 

In conclusion, in this paper, we have studied the twist-three fragmentation function 
contribution to the inclusive hadron's SSA in pp scattering p^p — > hX. With a simple 
parametrization for the twist-three fragmentation function, we estimate its contribution to 
the SSAs of 7T° production at RHIC energy. We find that the contribution of the twist-three 
fragmentation function is comparable to that of the twist-three distribution function from 
the polarized nucleon. We comment on the possibility to use our approach to describe the 
large difference of the SSAs between the r\ and 7r° meson. We emphasize that one need to 
include both contributions from twist-three distribution and fragmentation functions into a 
global analysis, in order to better understand the single spin asymmetry for the inclusive 
hadron production. 

We thanks Jianwei Qiu and Yuji Koike for helpful discussions. This work was supported 
in part by the U.S. Department of Energy under contract DE-AC02-05CH11231. We are 
grateful to RIKEN, Brookhaven National Laboratory and the U.S. Department of Energy 
(contract number DE-AC02-98CH10886) for providing the facilities essential for the com- 
pletion of this work. 
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